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Abstract 



< 

*i . In this communication I consider the applications of several kinds of approx- 

r^ ' imations of real functions to the problem of verified computation (reliable com- 

puting) of the range of implicitly defined real function x n+ i = G(x\, • • • , x n ), 
where dependency F(xi, ■ ■ ■ ,x n +x) = is defined on some compact domain 
by a sufficiently smooth real function F(xi, ■ ■ ■ , x n+ i). Constructive version 
of Kolmogorov-Arnold and implicit function theorems, results about floating- 
point approximation, floating-point approximations which give lower-bound and 
upper-bound estimates of some real functions, and approximate algebraic com- 
l/"~) , putation are used for the purpose. The rigorous theory can be build on the base 

^^ ' of analysis on manifolds over floating points domains. In the text we demon- 

Cm , strate our approach on examples. 

o 
o 



Introduction 



,—, ' In H the notion of Approximate Algebraic Computation (AAC) is formulated 

and shortly discussed. The subject matter of this talk lies in the area between 
C^ the geometry of numbers and the analysis of real functions and approximate 

algebraic computation of the functions. More specifically I want to discuss ap- 
proximate algebraic computation aspects of relation between the geometry of 
^A , parametric minima of convex and distance functions and the analysis of func- 

5^ ' tions which determine critical determinants. Let 



\ax + (3y\ p + \jx + 5y\ p < c\ det(aS ~ f3j)\ p/2 , 

be a diophantine inequality defined for a given real p > 1; hear a, (3, 7, d are real 
numbers with a8 — (3j ^ 0. H. Minkowski in his monograph [ |l3| raise the ques- 
tion about minimum constant c such that the inequality has integer solution 
other than origin. This Minkowski's problem can be reformulated as a conjec- 
ture concerning the critical determinant of the region | x \ p + \ y \ p < 1, p > 1. 
Mentioned mathematical problems are closely connected with Diophantine Ap- 
proximation. 



Also I want to proposed algorithms for approximate algebraic computation. 
By the computation of the algorithms A. Malishev and I have investigated 
the Minkowski's conjecture and proposed the strengthen Minkowski's analytic 
(MAS) conjecture. For verified computation (reliable computing) we used in- 
terval analyses. Methods and algorithms for interval evaluation of explicitly 
and implicitly defined real functions are used. In the paper we briefly consider 
following topics: 

critical determinant of a body and diophantine approximation; 
the problem; 

Minkowski's analytic conjecture; 
Kolmogorov- Arnold's theorem; 

computer experiments and strengthen Minkowski's analytic conjecture; 
interval-analytic methods; 
algorithms for approximate algebraic computation and implementation. 



2 Critical Determinant of a Body and Diophantine Approximation 

Critical determinant is one of the main notion of the Geometry of Num- 
bers. Recall the definitions [Q. Let V be a set and A be a lattice with base 
{a\, . . . ,a n } in R". A lattice A is admissible for body T> (T> — admissible) if 
V P| A = or 0. Let d(A) be the determinant of A. The infimum A(V) of deter- 
minants of all lattices admissible for T> is called the critical determinant of T>; 
if there is no I?— admissible lattices then puts A (I?) = oo. A lattice A is critical 
if d(A) = A(£>). 

Critical determinant is closely connected with diophantine approximation, solv- 
ing inequalities F(x\ , . . . , x n ) < c in integer numbers Xi,...,X n (with some 
restrictions, for instance, x — (xi, . . . ,x n ) =/= 0). Usually in the geometry of 
numbers the function F(x) is a distance function. A real function F(x) defined 
on R™ is distance function if 
(i) F{x) >0 ; xeR",F(0) = 0; 
(ii) F{x) is continuous; 

(iii) F(x) is homogenous: F(\x) = XF(x), A <E R. 

The problem of solving of diophantine inequality F(x) < c, with a distance 
function F are investigated. 

Let M be the closure of a set M and #P be the number of elements of a finite 
set P. An open set S C R n is a star body if S includes the origin of R™ and 
for any ray r beginning in the origin #(r n (M \ M)) < 1. If F(x) is a distance 
function then the set 

M F = {x : F(x) < 1} 

is a star body. 

One of the main particular case of a distance function is the case of convex sym- 
metrical function F(x) which with conditions (i) - (iii) satisfies the additional 
conditions 



(iv)F(x + y)<F(x)+F(y); 
(v) F{-x) =F(x). 



3 The Problem 



In considering the question of the minimum value taken by the expression 
\x\ p + \y\ p , with p > 1, at points, other that the origin, of a lattice A of determi- 
nant d(A), Minkowski |13|] shows that the problem of determining the maximum 
value of the minimum for different lattices may be reduced to that of finding 
the minimum possible area of a parallelogram with one vertex at the origin and 
the three remaining vertices on the curve \x\ p + \y\ p = 1. The problem with 
p = 1,2 and oo is trivial: in these cases the minimum areas are 1/2, v3/2 and 
1 respectively. Let D p C R 2 = (x, y), p > 1 be the 2-dimension region: 

\x\ p + \y\ p <l. 

Let A (Dp) be the critical determinant of the region. Using analytic parameter- 
ization Cohn H gives analytic formulation of Minkowski's conjecture. Let 

A(p,a) = {t + a)(l + rPyi {1 + aP)-p , (1) 

be the function defined in the domain 

D p : oo>p>l, 1 < a < a p = (2 P - 1)*, 

of the {p, a} plane, where a is some real parameter; here r = r(p, a) is the 
function uniquely determined by the conditions 

A p + B p = 1, < r < t p , 

where 

A = A(p, cr) = (1 + T p y L p - (1 + a p yp 

B = B(p,a) = <j(l + cr p yp +t(1 + t p )-p, 
t p is defined by the equation 

2(l-T p ) P = l+T p P , 0<Tp<l. 

In this case needs to extend the notion of parameter variety to parameter man- 
ifold. The function A(p, <r) in region D p determines the parameter manifold. 
Minkowski's analytic conjecture: 

For any real p and r with conditions p > 1, p =/= 2, 0<t<t p 

A(p,<r) > Tntn(A(p,l),A(p,ff p )). 



For investigation of properties of function A(p, a) which are need for proof 
of Minkowski's conjecture |OI, 0] we considered the value of A = A(p, a) and 
its derivatives A CT , A a2 , A p , A ap , A o on some subdomains of the domain 
Dp [S 10, O, O]. The analytical computation of the derivatives is a problem 
of computer algebra. 

4 The Theorem of Kolmogorov- Arnold 

Let F(x%, • • • , x n+ i) = 0, n > 1 be a sufficiently smooth real function of two 
variables. The rough form of the theorem of Kolmogorov- Arnold states that: 

Theorem (Kolmogorov, Arnold) 
Any sufficiently smooth real function can be represented as a superposition of 
functions of two variables. 

Let us demonstrate the constructive version of the theorem on examples of 
function A = A(p,a,r) and it's derivatives. At first expressing A a , A a2 , A p , A , A , 
in terms of a sum of derivatives of "atoms" s,; = o v ~ % , ti = t p ~ z , a,; = 
(1 + a*)"*"', h = (1 + rf)- l -p, A = b - a , B = rb Q + aa , a t = 
A p ~ l , (3i — B p ~ l (i — 0, 1, 2, . . .). Then by the implicit function theorem com- 
puting r = r(p, a) by means of the following iteration process: 

Ti+1 = (1 + 7?)* ((1 - ((1 + rf)-~v - (1 + a p )-*) P )* - <r(l + v p r*), 

For approximate computation of the expression for r p we apply the following 
iteration: 



(r P ) i+1 = 1 - (2-5)(l + (t^)5, p > l, ( Tp ) G [0,0.36]. 

So we really have represented the function A as the function A(p, a) of 
two variables. The same fact is true for it's derivatives Now we can compute 
expressions for A, A a , A ct2 , A p , A ap , Ao by means of approximate 
algebraic computations. 

5 Strengthen Minkowski's analytic conjecture 

Based on some theoretical evidences and results of mentioned computation 
A.V. Malishev and author proposed 

Strengthen Minkowski's analytic (MAS) conjecture: 



For given p > 1 and increasing a from to o~ p the function A(p, a) 

1) increase strictly monotonous if 1 < p < 2 and p > p^ 1 ', 

2) decrease strictly monotonous if 2 < p < p^ 2 ', 

3) has a unique maximum on the segment (1,<7 P ); until the maximum A(p,o~) 
increase strictly monotonous and then decrease strictly monotonous if p^ 2 ' < 
p < p (1) ; 

4) constant, if p = 2, 

here p^ 1 ' > 2 is a root of equation A ■ 2 \ <J=ry = 0; p^ 2 ' > 2 is a root of equation 

A^ 2 | ff=1 =0. 

It is seems that conjecture (MAS) did not proven for any parameter p except 

trivial p = 2. 



6 Interval-analytic methods 

Let X = (xi, • • • , x n ) = ([x 1; xi], • • • , [x n , x n ] be the n-dimensional real inter- 
val vector with Xj < Xi < Xj ("rectangle" or "box"). The interval evaluation of 
a function G(xi, • • • , x n ) on an interval X is the interval [G, G] such that for any 
x € X, G(x) € [G, G]. The interval evaluation is called optimal if G — minG, 
and G = max G on the interval X. 

In the communication I consider the case n = 2 R It is sufficient for Minkowski's 
conjecture. For the purpose we used modified variant of the method @ which 
we called 

Malyshev's Method: 

Let D be a subdomain of D p . Under evaluation in I? a mentioned function the 
domain is covered by rectangles of the form 

[p,p; s.,o\- 

In the case of formula (1) expressing A CT , A a2 , A p , A ap , Ao in terms of 
a sum of derivatives of "atoms" S{ = o- p ~ l , t% = t p ~ j , a* = (1 + a p )~ l ~~p , bi = 
(1 + rPy'-i, A = b - a , B = rb + aa Q: a, = A^\ fa = B*-* (* = 
0,1,2,...). one applies the rational interval evaluation to construct formulas 
for lower bounds and upper bounds of the functions, which in the end can be 
expressed in terms of p, p, a, a, r, r, ; here the bounds r, r, are obtained 
with the help of the iteration process: 



(1 +£)*((! -((l+#)-5-(l + a£) *~) l r~ -a(l + o£) £), 



t s ; +1 = (l+f)I((l-((l+f)^-(l+^)-?F)5-a(l + aP)-f). 



x By the result of A. Kolmogorov and V. Arnold, any sufficiently smooth real function can 
be represented as a superposition of functions of two variables. 



i = 0, 1, • • • 

As interval computation is the enclosure method, we have to put: 

[Z, t] = [t N , t N ]f][T , To] . 

N is computed on the last step of the iteration. 

For initial values we may take : [t 0l to] — [t , to] = [0, 0.36]. 

7 Algorithms for approximate algebraic computation 

Below each of 6 first algorithms has 4 different forms: 
(i) approximate algebraic computation of the given expression in a given point 
at the floating point representation; 

(ii) approximate algebraic computation of a lower-bound estimate of the given 
expression in a given point at the floating point representation; 
(iii) approximate algebraic computation of a upper-bound estimate of the given 
expression in a given point at the floating point representation; 
(iv) approximate algebraic computation of an interval evaluation of the given 
expression over given intervals at the floating point representation; 
Here we give names, input and output of algorithms for interval evaluation only. 
But all algorithms for Minkowski conjecture are implemented and tested. 

Algorithm TPV 

Input: An implicitly defined function t p from Paragraphs 3 and 6. 

\p,p; 3.M- 

Method: Shortly described in Paragraphs 3 and 6. 
Output: The interval evaluation of t p . 

Algorithm TAUV 
Input: Implicitly defined function t from Paragraphs 3 and 6. 

[p,p; 3.M- 

Method: Shortly described in Paragraphs 3 and 6. 
Output: The interval evaluation of t. 

Algorithm LOV 
Input: Function 1° — A(p, a) — A p from Paragraphs 3 and 6. 

[p,p; g_,a}. 

Method: Shortly described in Paragraphs 3 and 6. 
Output: The interval evaluation of 1°. 

Algorithm LI V 
Input: Function I 1 = A(p, a) — A p from Paragraphs 3 and 6. 
[P,F, Qi,a]. 
Method: Shortly described in Paragraphs 3 and 6. 



Output: The interval evaluation of I 1 . 

Algorithm GV 

Input: A function g(p, er) which has the same sign as function A CT . 

[p,p; Qi,W). 

Method: Shortly described in Paragraphs 3 and 6. 

Output: The interval evaluation of g(p,a). 

Algorithm HV 

Input: A function h{p, a) which is the partial derivative by a the function 

fffecr). 

[P,P; 2., a]. 

Method: Shortly described in Paragraphs 3 and 6. 

Output: The interval evaluation of h(p, a). 

Next two algorithms are described in M. 

Algorithm MonotoneFunction 
Input: A real function F(x, y) monotonous by x and by y. 
Interval \x, x; y, y] . 
Output: The interval evaluation of F. 

Algorithm RationalFunction 
Input: A rational function R(x,y). Interval \x,x;y,y\. 
Output: The interval evaluation of R. 



Implementation 

Algorithms are implemented on PL/1(0) and (partly) on C and C++. There 
is also an implementation of some of the algorithms on Reduce3.3-3.X in T. 
Sasaki's package Q of arbitrary precision rational arithmetic. 
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